A note on the extension of the Dirac method 
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In this note we extend the Dirac method to partial differential equations involving higher order 
roots of differential operators. 
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In deriving his celebrated equation, the Dirac stroke of genius was to use the anti-commuting nature of the Pauli 
matrices 

{a^&k} = 2S jk , ft* = 1,2, 3) (1.1) 

to rewrite the Pythagorean theorem as follows 



A = VB 2 +C 2 =B& 3 +Ca k (j^k) (1.2) 

where A,B,C can be either numbers or operators. Eg II .21 allows to eliminate the square root from the relativistic 
Schrodinger equation 



>• ! ihd t ^ = c \/p 2 + m 2 c 2 ^ , (1.3) 

00 

C ■ 1 which, can accordingly be written as 
00 ' 

■ ihd t £ = c(p&i + mca 3 ) * . (1.4) 

where is a 2-component vector. This equation resembles the Dirac equation, but it should not be confused with it 
or with the Pauli equation. It does not involve, indeed, 4-component spinors and 4x4 matrices, and is interpreted 
as the equation ruling the free evolution of a relativistic spinless particle. The underlying physical problems have 
already been discussed in Refs. |2j, l3j, and will not be reconsidered here. 

What we are concerned is whether it is possible to generalize eq. (11.21) as follows 



A = t/B 3 + C 3 = Bfj+C f k (j ^ k) (1.5) 

with the matrices fj (j — 1,2,3) satisfying the conditions 

ff = 1 t 2 r k + tj h fj + f k fj = . (1.6) 

If do exist, these matrices have to be associated with the cubic roots of unity 
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that satisfy the conditions 

e 3 = l, e + e+ + e-=0, e+ e_ = 1 
We construct two of the possible matrices as 



1 
fi = I 1 



e + 
t-2 = | e_ 

1 



1 

which are easily shown to satisfy the conditions (|1.6[) . Accordingly the equation 



dt$= ZJdl + k $ 

with k being a constant, can be cast in the form 

= (d x h +fcf 2 )t 



t(* = 0) = lo , 



(1.8) 



(1.9) 



(1.10) 



(1.11) 



where $ is a 3-component vector. This increase of dimensionality is the price to be paid if we want to use a Dirac-like 
procedure to avoid the use of operators raised to a non integer number (further comments on this point can be found 
inRef. 0). ^ 

The solution of eq. (jl.lip can be written as 



$ = cxp {t (d x h+k f 2 )} $ . 



(1.12) 



The exponential operator in the last equation cannot be naively disentangled because the matrices f are not com- 
muting. Indeed, one has 



[Tl,T 2 \ 



T3 



with 



T3 



1 

e+ 
e_ 



(1.13) 



(1.14) 



Also this matrix satisfies the conditions (|1.6[) . 

The explicit solution of our problem can be found using operator ordering methods 4]. However, just to provide a 
feeling on the possibilities offered by the present method, we use the Trotter formula [5j 



exp 



(A + b) 



= lim 



exp 



n 



exp 



n 



to write cq. (|1.12[) as follows 



$ = lim 

n— >oo 



tk 

exp [ — t% J exp 



td x 



(1.15) 



(1.16) 



By considering small time intervals 5 = t/n, and evaluating the effect of the exponentials step by step, we get an 
iterative solution in the form 



<J> n = exp (<5 k f 2 ) exp (6 d x n) 



(1.17) 



The action of the evolution operator on the function on its right can be easily obtained by noting that the exponentials 
can be expanded in finite terms, as a consequence of the cyclical properties of the matrices f, namely f 3m = 1, thus 
getting @ 



exp (a Tj) = C (a) 1 + C\ (a) fj + C 2 (a) f ? 



with 



00 n 3l+m 



1=0 



(31 + to)! 



(to = 0,1,2). 



(1.18) 



(1.19) 
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These functions are recognized as the pseudo- hyperbolic functions introduced in Q, or, in more general terms, as the 
generalized trigonometric functions discussed recently in Refs. @, 

According to the previous relations, at the first iteration we get ($ = (<pi(x),cf2(x), ipz{x)) T ) 

= cxp (S k f 2 ) exp (S d x n) J> 




C (6k) e+diSk) C 2 (6k) \ fC (8d x )ip 1 {x) + C 1 (5d x )ip2(x) + C 2 {8d x )<p 3 (x)\ 
e„C 2 (5k) C {Sk) e-d(8k) C 2 (5 d x ) ^(x) + C (5 d x ) (p 2 (x) + C^S d x ) <p 3 (x) . (1.20) 
diSk) e+C 2 (Sk) Co(Sk) J Vd(5a x ) < Pl (x)+ C 2 (5d x )ip 2 (x) +C (6d x ) <p 3 (x) J 



The successive iterations are straightforwardly obtained along the same lines. 

It is evident that the method we have outlined open new perspectives and is amenable for further implementations 
which will carefully be discussed in a forthcoming investigation where the case of higher order roots of unity will be 
discussed. 



[1] P. A. M. Dirac, Proc. R. Soc. London A 117, 610 (1928); ibid. 118, 351 (1928). 

[2] Kunling Wang, Tao Liu, Mang Feng, Wanli Yang, Kelin Wang, Phys. Rev. A 82, 064501 (2010) 

[3] D. Babusci, G. Dattoli, M. Quattromini, arXiv:1101.5066vl[math-ph]. 

[4] G. Dattoli, P. L. Ottaviani, A. Torre, and L. Vazquez, Riv. Nuovo Cimento 20, n. 2 (1997). 

[5] H. F. Trotter, P. Am. Math. Soc. 10, 545 (1959). 

[6] G. Dattoli, M. Migliorati, P. E. Ricci, arXiv: 1010 . 1676vl [math-ph] . 

[7] P. E. Ricci, Le Funzioni Pseudo-Iperboliche e Pseudo-Trigonometriche, Pubbl. Istit. Mat. Appl. Fac. Ing. Univ. Stud. Roma 
12, 37 (1978). 

[8] R. M. Yamaleev, Adv. Appl. Clifford Al. 10 (1), 15 (2005). 
[9] D. Babusci, G. Dattoli, E. Di Palma, E. Sabia, arXiv: 1103. 2321vl [math. CA] . 
[10] K. Fujii, H. Oike, [arXiv:quant-ph/0604115vl]. 



